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Soil salinity problems are widespread around the globe with increased risk of spreading over the years. The 
fungus Piriformospora indica, identified in Indian Thar desert, colonizes the roots of monocotyledon plants 
and provides resistance towards biotic as well as abiotic stress conditions. We have identified a cyclophilin 
A-like protein from P. indica (PiCypA), which shows higher expression levels during salinity stress. The 
transgenic tobacco plants overexpressing PiCypA develop osmotic tolerance and exhibit normal growth 
under osmotic stress conditions. The crystal structure and NMR spectroscopy of PiCypA show a canonical 
cyclophilin like fold exhibiting a novel RNA binding activity. The RNA binding activity of the protein and 
identification of the key residues involved in the RNA recognition is unique for this class of protein. Here, we 
demonstrate for the first time a direct evidence of countering osmotic stress tolerance in plant by genetic 
modification using a P. indica gene. 

Abiotic stress leads to a series of morphological, physiological, biochemical and molecular changes that 
adversely affect plant growth and productivity 1,2 . High salinity is one of the major threats to plant growth 
and productivity; it causes both hyperionic and hyperosmotic stress and can lead to plant death 3 . Plants 
basically counteract these negative effects of salinity stress by activation of various biochemical and metabolic 
pathways 4 that include accumulation of osmolytes like proline, fructose, myo-inositol, betaine, sucrose etc. 5 , 
maintaining the intracellular ion homeostasis and scavenging of reactive oxygen species (ROS) 68 . The over- 
expression of various plant genes has been shown to provide salinity stress tolerance to plants 3 ' 9,10 . 

P. indica, a plant-root-colonizing basidiomycete fungus, was discovered in the Indian Thar desert and has been 
shown to provide strong growth-promoting activity during its symbiosis with a broad spectrum of plants 11 . 
Symbiotic relationship of P. indica with plants results in better plant growth, higher seed yield, early flowering, 
and biotic and abiotic stress tolerance responses 1213 . It has been also reported that mutual nutrients like phosphate 
transport take place by P. indica 1 *. Salt-stress studies have shown positive effects of P. indica on barley 12 . P. indica 
also provides resistance towards fungal diseases and tolerance to salt stress in monocotyledonous plants 15 . 
However, genes from this beneficial fungus, P. indica, in providing salinity stress tolerance to plants have not 
been exploited yet. 

Cyclophilins (Gyps) are highly conserved, ubiquitous in nature and catalyze the interconversion of peptidyl 
prolyl imide bonds in peptide and protein substrates 16 . In eukaryotic cells, Cyps are found in all cellular compart- 
ments with a variety of functions being ascribed to them including cell division, transcriptional regulation, protein 
trafficking, cell signaling, pre-mRNA splicing, molecular chaperone mechanism and stress tolerance 16 . A striking 
feature of plant Cyp genes is that their expression is induced in response to various stresses as well as abscisic acid 
(ABA). Further, wide distribution and ubiquitous nature of Cyps signifies their fundamental importance in plant 
survival. Although diverse functions of Cyps have been suggested in plants, their physiological relevance and 
the molecular basis of stress-responsive expression are still largely unexplored 17 20 . The role of plant Cyp in 
conferring salt-stress tolerance in rice has been reported 21 . Recently, we have identified 28 rice, 35 Arabidopsis 
and 8 yeast Cyps genes from their respective genomes through a genome-wide analysis 16 . Earlier, we have reported 
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sequence- specific 'H, 13 C and 15 N resonance assignments for a 
Cyclophilin A-like protein from Piriformospora indica (PiCypA) 
which was found to be upregulated during osmotic stress conditions 22 . 

Though Cyps are known to play role in pre-mRNA splicing, their 
RNA-binding activity has not been well described. RNA-binding 
proteins (RBPs) play major role in RNA metabolism, pre-mRNA 
processing, transport, stability and translation 23,24 . RBPs are involved 
in posttranscriptional regulatory machinery that further allows 
plants to respond against environmental stresses 25 . Many RBPs have 
been reported to be regulated during abiotic stresses in plants like 
MA16, a glycine rich RBP which is induced by dehydration stress and 
ABA in maize embryos 26 , Tudor-SN (TSN) protein which is highly 
expressed during drought stress in Arabidopsis 27 . DEAD-box RNA 
helicase (ZmDRHl) in Zea mays 2S and glycine rich RBP from 
Sorghum bicolor 29 have been reported to be upregukated under abi- 
otic stress. Here we report the crystal structure of PiCypA protein, 
the role of PiCypA in imparting osmotic stress tolerance to plant and 
the unique property of PiCypA to bind RNA. Some of the interacting 
partners of PiCypA protein were also upregulated under salt stress in 
the transgenic plants. Furthermore this tolerance is observed in T 2 
generation of transgenic plants and the possible mechanism of the 
tolerance is by maintaining the antioxidant machinery of plants 
under stress conditions. 



Results 

Identification of PiCypA and its overexpression in tobacco plants. 

P. indica was grown at different salt concentrations and good growth 
of P. indica was observed at 200 and 400 mM NaCl, while there was 
moderate growth at 600 mM NaCl and no growth at 800 mM NaCl. 
To identify the genes involved in the high salt tolerance, the 
functional screening of salinity tolerant genes was performed. For 
this purpose E. coli cells were transformed with phagemids derived 
from a cDNA library (made from RNA isolated from P. indica grown 
at 600 mM NaCl) as described in materials and methods. A total of 
over one million recombinant bacterial cells (E. coli SOLR cells) were 
plated on LB agar medium supplemented with 600 mM NaCl and 
without salt as control. A lawn of colonies on the control and only 36 
colonies on the plate containing 600 mM NaCl were observed. Out 
of 36 salt tolerant colonies we selected one for further analysis, 
which was found to be a gene encoding cyclophylin A-like protein 
(PiCypA; accession number GQ214003). Functional significance 
of the PiCypA gene was established by overexpressing the com- 
plete open reading frame (ORF; 495 bp) in tobacco plants using 
Agrobacterium-medmted transformation with the gene construct 
pCAMBIA1301-PtC>pA (Fig. la). 

The same vector not containing the PiCypA gene called empty 
vector control (VC), pCAMBIA1301, was also transformed to plants 
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Figure 1 | Analysis of PiCypA overexpressing transgenic tobacco plants (T x generation) under osmotic stress, (a) Structure of T-DNA region of 
pCAMBIA1301 containing PiCypA gene (495 bp) in Xhol and EcoRI restriction enzyme sites of MCS region with CaMV35S promoter and poly A 
terminator (pCAMBIA1301-PiCypA). (b) Visualization of fi-glucuronidase (GUS) activity in the leaf tissue of WT and vector control (VC) controls and 
the transgenic lines (L12.1, L12.2 and L12.3) by using X-gluc. The hygromycin resistant putative tobacco transformants were used for GUS assay, (c) PCR 
analysis of the transgenic lines by the use of PiCypA gene specific forward and reverse primers. Lane 1 (M) is DNA marker, lane 2 is WT, lane 3 is VC and 
rests of lanes (4—6) are transgenic lines showing the required amplification (495 bp), (d) The PiCypA transgenic lines and control plants (WT and VC) 
were analyzed by the Southern blot hybridization to show the integration of the PiCypA gene along with CaMV35S promoter and poly A terminator 
(1.1 kb) in the transgenic lines. The genomic DNAs were digested with Hindlll and probed with radiolabeled ORF of PiCypA gene. The gel is cropped for 
clarity, (e) Relative mRNA expression in all the three transgenic lines checked by Real Time PCR. (f) Western blot analysis of the WT and VC control 
plants and T] transgenic lines (12.1, 12.2 and 12.3) using PiCypA polyclonal antibodies. The gel is cropped for clarity. Full gel figures for (d) and (f) are 
shown in supplementary figure S4. 
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treated as an additional control along with the wild type ( WT) plants. 
Seeds from the T 0 plants, when plated onto hygromycin-containing 
medium, segregated in the expected 3 : 1 ratio of HygVHyg" (Table 1). 

Analysis for the PiCypA transgene integration was carried out on 
the Tj generation of transgenic plants. We selected 16 independently 
transformed T x transgenic plants and grew them to maturity; how- 
ever, only 3 different lines (L12.1, L12.2 and L12.3) were randomly 
selected for subsequent analysis. The hygromycin resistant putative 
tobacco transformants and VC plants were confirmed by /i-glucur- 
onidase (GUS) assay and found to be positive whereas WT was GUS 
negative (Fig. lb). Integration of PiCypA gene (495 bp) was con- 
firmed in the GUS positive transgenic lines by PCR using gene spe- 
cific primers, whereas WT and VC plants were PCR negative 
(Fig. lc). The morphological and growth characteristics of T l plants 
were similar to the WT and VC control plants. The stability of 
integration of the transgene was further confirmed by probing of 
Hindlll digested genomic DNA with an equal amount of PiCypA 
cDNA. The expected band pattern per insertion was observed at 
1.1 kb position (containing the entire cauliflower mosaic virus pro- 
moter and full ORF of PiCypA cDNA) (Fig. Id, lanes 2-4). No band 
was detected in the WT plants (Fig. Id, lane 1). Real time PCR of all 
the three lines showed the expression of transcript in all of the trans- 
genic lines (Fig. le). Western blot of the Ti transgenics by PiCypA 
polyclonal antibodies identified an 18 kDa band corresponding to 
the purified recombinant protein, and its intensity varied between 
the different transgenic lines (Fig. If, lanes 3-5). A very faint protein 
band of the same size was also observed in the WT and VC plants 
(Fig. If, lanes 1 and 2). This could be possibly due to the cross- 
reactivity of the polyclonal antibodies with native tobacco protein. 
All these findings confirmed the stable integration of the PiCypA 
gene into the transgenic plants. 

PiCypA provides tolerance to excess osmotic stress. To test for 
osmotic stress tolerance, leaf disks from T] transgenic, control 
tobacco plants (WT and VC) were floated separately on 200 mM 
NaCl and H 2 0 for 72 h. Salinity-induced loss of chlorophyll was 
lower in PiCypA overexpressing lines (L12.1, L12.2 and L12.3) 
compared with those from the WT and VC plants (Fig. 2a). The 
damage caused by stress was reflected in the degree of bleaching 
observed in the leaf tissue after 72 h. It was evident that transgenic 
lines had a better ability to tolerate osmotic stress. The measurement 
of the chlorophyll content of the leaf disks from different transgenic 
lines stressed with 200 mM NaCl provided further support for a 
positive relationship between the overexpression of PiCypA and 
tolerance towards osmotic stress (Fig. 2b). To assess the tolerance, 
we germinated the seeds of all three transgenic lines along with WT 
and VC in 200 mM NaCl and in water as control. All the above seeds 
germinated well in water (Fig. 2c, first row), while a better and high 
germination rate of transgenic lines (12.1, 12.2 and 12.3) was 



observed as compared to WT and VC (Fig. 2c, bottom row). The 
T t seedlings were morphologically similar when grown with or 
without NaCl. To assess the effect of high salt (200 mM NaCl) on 
the growth and development of PiCypA overexpressing lines, WT 
and VC plants, the three week old seedlings (after growing in water) 
were further grown in the soil continuously irrigated by 200 mM 
NaCl. In the presence of NaCl, the WT and VC plants showed 
growth retardation, whereas transgenic plants did not develop any 
sign of stress (Fig. 2d, bottom row). It was observed that without salt 
all plants had similar growth (Fig. 2d, upper row). The T l transgenic 
plants grew normally and set viable seeds under continuous osmotic 
stress, whereas the WT and VC plants could not survive till maturity 
under continuous osmotic stress (Fig. 2e). Overall, the results 
showed that the PiCypA overexpressing Ti transgenic lines have 
shown better ability to tolerate osmotic stress. This is the first 
direct evidence to show osmotic stress tolerance in plant by a P. 
indica fungus gene. 

Growth parameters and viable seeds under osmotic stress in 
PiCypA-transgenic plants. The percent seedlings survival of trans- 
genic lines grown in 200 mM NaCl did not significantly change as 
compared to water grown WT and VC controls (Table 1). Several 
critical growth parameters, such as plant height, fresh weight of 
leaves, total protein content, time required for flowering and seed 
weight were scored as an indicator of osmotic stress to test tolerance 
in Tj transgenic lines (Table 1) as compared with the WT plants 
grown in H 2 0. It should be noted here that similar data for WT/VC 
plants grown under osmotic stress could not be obtained because 
these plants failed to sustain growth in the presence of 200 mM NaCl 
(Fig. 2e). Overall performance or total seed yield of the transgenic 
plants was unaffected in the presence of 200 mM NaCl. The 
flowering time in osmotic-stressed plants was delayed by a week 
when compared with plants grown without stress, but there was no 
change in the flower morphology and yield (Table 1). Significantly, 
the pod and seed size, seed number and seed weight of salt grown 
transgenic lines were statistically similar to those from WT and VC 
plants grown without osmotic stress (Table 1). 

PiCypA overexpressing tobacco lines confers tolerance to oxida- 
tive stress. In response to 200 mM NaCl, PiCypA transgenic lines 
exhibited lesser accumulation of malondialdehyde (MDA) (lipid 
peroxidation product), H 2 0 2 and ion leakage as compared to 
controls seedlings (WT and VC), suggesting more tolerance to oxida- 
tive stress in transgenic tobacco under the osmotic stress (Supple- 
mentary fig. SI a-c). We observed higher increase in the activities of 
Catalase (CAT), ascorbate peroxidases (APX) and glutathione 
reductase (GR) in transgenic lines as compared to control plants 
(WT and VC) under osmotic stress (Supplementary fig. SI d, e 
and f). Moreover, the proline accumulation is another common 



Table 1 | Comparison of segregation ratio, plant seedlings survival (%) and various growth parameters of the WT, VC controls and 
Cyclophilin T] transgenic plants grown in the presence of water and 200 mM NaCl, respectively 

Parameters Water grown WT** Water grown VC Salt grown LI 2.1 Salt grown LI 2.3 Salt grown LI 2.7 



Segregation ratio, Hr:Hs [n]* 


0 


3.2: 1 (118) 


3.12: 1 (158) 


2.9: 1 (143) 


3.18: 1 (133) 


Plant seedling survival (%) 


96 ± 1 


95 ± 2 


92 ± 2 


97 ± 2 


96 ± 2 


Total protein mg/g FW 


0.78 


0.76 


0.79 


0.83 


0.80 


Fresh weight/ g/leaf 
Height/ cm 


2.8 


2.9 


2.91 


3.01 


3.2 


107 


110 


1 10 


112 


109 


Flowering time, days 


1 16 ± 2 


122 ± 3 


102 ± 3 


1 13 ± 2 


108 ± 2 


Seed weight per Pod (grams)/ 


140 


132 


142 


158 


147 



'Recordings made from seed. Experiments were repeated at least three times, and the values represent means of different samples. SD in each case was <3%. In all cases, n = 1 5. n is mentioned for each line 
in brackets. 

"WT plants did not survive under osmotic stress. 
Recordings made from mature leaf. FW, fresh weight. 
Recordings made from seed. 
^Recordings made from mature plant. 
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Figure 2 | Analysis of osmotic stress tolerance in the transgenic lines, (a) Leaf disk senescence assay for osmotic tolerance in transgenic tobacco plants 
(L12.1, L12.2 and L12.3) in 200 mM NaCl. Transgenic lines showing high osmotic stress tolerance as compared to controls (WT and VC). (b) Remaining 
chlorophyll content from leaf disks of WT, VC plants and the transgenic lines after incubation 200 mM NaCl solution (Clh A and Clh B are shown in dark 
red and light green bar respectively). Disks floated in water served as the experimental control (yellow bar represents total chlorophyll content in the 
control condition), (c) Germination and seedling growth ofWT, VC and transgenic lines (L12.1, L12.2 and L12.3) under water condition (for control) 
and 200 mM NaCl (for experiment), (d) Phenotypic comparison of growth of T] transgenic tobacco plants overexpressing PiCypA gene (L12.1, L12.2 
and L12.3) with controls plants (WT and VC) in the presence of water and 200 mM NaCl. Seeds were germinated for 7 days and plantlets grown on 
water (upper row) and on 200 mMNaCl (lower row), (e) WT, VC and transgenic plants (L12.1,L12.2 andL12.3) in soil pots supplied with 200 mMNaCl 
solution. Note that the WT and VC plants could not sustain growth under osmotic stress, (f) Phenotypic comparison of growth of T2 transgenic tobacco 
plants overexpressing PiCypA gene (L12.1.1, L12.2.1 and L12.3.1) with control plants (WT and VC) in the presence of water and 200 mM NaCl. 



adaptive mechanism for providing osmotolerance to plants against 
abiotic stress 42 . In this study we observed that the proline accumu- 
lation and relative water content were also maintained in transgenic 
lines as compared to control plants (WT and VC) (Supplementary 
fig. SI g and h). 

T 2 transgenic tobacco plants grow normally and set viable seeds 
under osmotic stress. The hygromycin resistant Tj transgenic seeds 
were grown in soil under normal conditions. These plants (T 2 
generation) were morphologically similar to that of WT and tested 
for osmotic stress tolerance. Osmotic stress-induced loss of chloro- 
phyll in leaf disk assay was lower in PiCypA overexpressing T 2 lines 



(L12.1.1.L12.2.1 and L 1 2 . 3 . 1 ) compared with those from the WT and 
VC plants (Supplementary fig. S2 a and b). PiCypA overexpressing 
T 2 lines under osmotic (200 mM NaCl) stress accumulated less H 2 0 2 
in comparison with WT and VC as tested by DAB staining method 
and its quantitation (Supplementary fig. S2 c and d). The 
concentration of proline, an important osmoregulator, was higher 
in PiCypA overexpressing T 2 lines as compared to WT and VC 
(Supplementary fig. S2e). The growth retardation was also re- 
ported in three weeks old WT and VC plants whereas T 2 trans- 
genic lines grew happily under continuous salt irrigation (Fig. 2f, 
lower row). It was observed that without salt all the plants showed 
similar growth (Fig. 2f, upper row). 
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Molecular architecture and dynamics of PiCypA. To understand 
the molecular basis of PiCypA's role in providing salt stress tolerance, 
we performed structural and dynamics studies. Crystal structure of 
PiCypA was solved at 1.97 A resolutions by molecular replacement 
method. Crystallization and data collection statistics have been 
described previously 30 . The model includes the entire polypeptide 
chain and was refined to R wor k an d Rf ree values of 0.159 and 0.183, 
respectively (Supplementary Table SI). Crystal arrangement has 
shown three molecules per asymmetric unit arranged in a side-by- 
side orientation. PiCypA is a monomeric single domain protein and 
forms a canonical cyclophilin fold comprising two a-helices and 
eight /i-strands (Fig. 3a). Eight /i-strands form a flattened cylinder 
with /i-barrel like fold where two helices sit on top and bottom of the 
barrel. Six phenylalanines (Fll, F25, F39, F56, F115 and F132) and 
one tyrosine (Y51) arranged in edge-to-face orientation form a tight 
hydrophobic core. Around this core, other hydrophobic aliphatic 
residues (V9, 113, 123, L27, L42, 159, L65, L101, M103, 1117 and 
V135) form a closed cage structure (Fig. 3b). Eighty five residues 



comprising -52% of a full polypeptide chain (N45-F56, G68-P98, 
A104-Q114 and T119-V131) are present in five loop regions that 
have similar B-factor compared to average B-factor of 29.5 A 2 for the 
entire polypeptide chain. The rigid nature of loops in solution was 
also confirmed from the NMR dynamics data. The residues forming 
these loops have similar order parameter like those forming a-helices 
and /?-strands that confirms that the loops are rigid. We calculated 
order parameter by analyzing NMR relaxation parameters R lt R 2 , 
and 15 N-{'H} nOe. Complete sequence-specific NMR resonance 
assignments of the protein have been reported previously 22 . We 
had observed that chemical shifts of amide protons of L34 
(10.67 ppm), D88 (11.36 ppm) and H95 (10.74 ppm) and amide 
nitrogen of E89 (132.8 ppm) are unusually downfield shifted. The 
unusual downfield shifts arise due to deshielding effect caused by 
delocalized re-electron ring currents of the proximal aromatics 
rings of residues Y82, F86, F91 and W124 (Fig. 3c). 

As PiCypA shares 62% sequence homology with human cyclophi- 
lin A, we investigated two characteristic features of cyclophilin-like 
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Figure 3 | Molecular architecture and RNA binding activity of PiCypA. (a) /J-barrel like fold containing eight /?-strands and two helices forming cap at 
either ends (b) aromatic (cyan) and aliphatic (magenta) forming hydrophobic core (c) Residues showing unusual chemical shifts (Magenta) in NMR 
spectra and aromatic residues (cyan) responsible for it. (d) Western blot of PiCypA protein probed with anti-his antibody (upper panel) as loading 
control. For this 1 ug of BSA and increasing concentration of PiCypA (1, 2 and 4 pg) proteins were dot-blotted on precharged PVDF membrane. An 
identical blot was probed with a 13 mer radio labeled RNA oligo (5'-AUAGCCUCAACCG-3') and developed, (e) Isothermal titration calorimetry 
indicating hyperbolic nature for PiCypA binding to the 13 mer RNA (f) Selected regions from the superimposed 2D [ 15 N, 'H] HSQC spectra of free and 
complexed PiCypA shown in red and blue, respectively. The resonance assignments of residues with large chemical shift changes upon RNA binding (1:1 
molar ratio) are shown with arrows, (g) Combined chemical shift changes of PiCypA upon RNA binding in 1 : 1 molar ratio. Red line indicates one 
standard deviation for the mean chemical shift change of the dataset. Combined chemical shift changes were calculated by taking the square root of the 
sum of the squares of the chemical shift changes for the backbone 'H N and 15 N resonances after scaling the 15 N shifts by 0.2 (h) Mapped RNA binding 
region (Rl, R2 and R3) on the crystal structure of PiCypA. 
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fold containing proteins, PPIase (peptidyl prolyl cis-trans isomerase) 
domain architecture and cyclosporin A (CsA) binding activity. 
PiCypA carries a PPIase domain with all conserved residues (R58, 
F63, M64, Q67, A104, Fl 16, W124, L125, and H129). NMR titrations 
with CsA have shown that PiCypA has same conserved CsA binding 
surface (H57, R58, 160, Q66, G67, G68, H95 and 1117). Along with 
these residues, three more residues that are mutated from the human 
counterpart (C55T, C64M and C118T) completes CsA interacting 
region (Supplementary fig. S3 a-c). 

To explore the effect of high salt on the structure of PiCypA we 
performed salt titration experiments. In the 2D [ 15 N, 1 H] HSQC 
spectra except for minor changes in chemical shifts of surface 
exposed residues T55, S80 and K85 no significant backbone amide 
chemical shift perturbations (CSPs) were observed confirming that 
structural integrity of the protein is maintained even at 600 mM 
NaCl 22 . 

PiCypA specifically binds to RNA. We observed a novel RNA 
binding property of PiCypA when the protein was spotted on a 
charged PVDF membrane in increasing concentration. PiCypA 
showed strong RNA binding as shown in autoradiograph. Bovine 
serum albumin (BSA) was used as a negative control (Fig. 3d). 

To investigate specificity of RNA binding, PiCypA was titrated 
against different RNA constructs and we obtained quantitative ther- 
modynamic parameters of binding. The selected RNA sequences 
were of different length and sequences (1. 5'-GCU GCG-3'; 2. 5'- 
CGU GCG-3'; 3. 5'-AUA GCC UCA ACC G-3'; 4. 5'-GCU GCU 
GCG-3'). No RNA binding was observed for 6 mer or 13 mer poly- A 
RNA sequence. A 13 mer RNA sequence (AUAGCCUCAACCG) 
was the strongest binder to the protein. The binding constant (K d ) 
of 22.5 uM indicates moderate to strong binding with AH and AS 
values of —118.4 Kcal/mol and —369 cal/mol/deg respectively. The 
observed isotherms indicated a moderate to strong nature of binding. 
Therefore, corresponding stoichiometry of the binding could not be 
determined as the observed 'c value' of isotherm was 0.2. As the 
observed 'c value' is outside of the desired range of 1 < c < 1000 a 
hyperbolic isotherm as compared to sigmoidal isotherm was 
observed that does not allow correct determination of stoichi- 
ometry 3132 (Fig. 3e). 

13 mer RNA (AUAGCCUCAACCG) was titrated against U- 15 N 
labeled PiCypA at pH 6.5 and shifts in the backbone 15 N and 'H N 
resonances in 2D [ 15 N, 'H] HSQC spectra were monitored. A gradual 
shift in the resonances was observed for all the residues that showed 
changes indicating a fast exchange (t cx < Lis) between free and RNA 
bound PiCypA (Fig. 3f). No significant changes in chemical shifts 
were observed after 1 : 1 molar ratio of protein and RNA indicating 
stoichiometry of 1:1 in the complex hence indicating one RNA 
molecule binding to one PiCypA molecule. The residues showing 
the high cumulative backbone CSPs are most likely to be involved in 
the RNA binding. The residues showing significant CSPs are clus- 
tered into three major regions of the polypeptide chain. These 
regions are R; (L65, G68, T71, G76, T77, G78, G79 and R80), R 2 
(M103, A104, N105, G107 and T110), and R 3 (W124, L125, D126, 
G127, K128, H129 and V130) along with two individual residues T55 
and F115 (Fig. 3g). We mapped these residues on the crystal struc- 
ture of PiCypA. These three regions are clustered around the groove 
that comprises of aromatic as well as putative hydrogen bonding 
residues, making it more likely to interact with nucleic acid bases 
(Fig. 3h). 

Interaction of PiCypA with other stress induced proteins and 
their expression level. We performed yeast-2-hybrid assays to find 
putative interacting protein partners from P. indica cDNA library. A 
total of 15 interacting proteins of PiCypA were obtained and their 
putative functions are mentioned in Supplementary Table S2. Out of 
these 3 interacting proteins were found to be abiotic stress induced as 
reported earlier 933 36 . The expression level of these genes namely 



phosphoglycerate kinase (PGK), heat shock protein 6 (HSP6) and 
eIP-4A were found to be higher in transgenic lines as compared to 
WT (Fig. 4a-c) suggesting that PiCypA might be playing vital role in 
maintenance of mRNA integrity. We also performed semi- 
quantitative PCR to check the expression level of PGK, HSP6 and 
eIF-4A mRNAs (Fig. 4d). The high expression level of these stress- 
induced genes together with PiCypA provides increased tolerance to 
osmotic stress. 

Discussion 

Soil salinity problems have hampered cultivation of many crop 
plants in large parts of the world. Unregulated large-scale irrigation 
programs and poor soil management are major reasons for its 
increased risk of spreading over the years. High salinity stress con- 
ditions affect the cellular homeostasis and cause disturbance in water 
potential and ion distribution. This disturbance in water and ion 
potential leads to molecular damage, growth inhibition and even 
death of plants. The molecular chaperones from plants such as heat 
shock proteins (HSPs) and cyclophilins have been shown to play role 
in abiotic stress tolerance 353 738 . Cyclophilins have endogenous con- 
served cytosolic peptidyl-prolyl isomerase (PPiase) activity, which 
catalyzes the inter-conversion of the cis and trans conformation of 
peptide bonds of proline residues in polypeptides. This inter-conver- 
sion is a rate-limiting step during protein folding 39,40 . Cyclophilins 
have been reported to play important role in different cellular func- 
tions such as inhibition of cell proliferation, induction of apoptosis 
and suppression of migratory/invasive capacity 41 . Interestingly, in 
silico analysis have predicted Cyclophilin to be an abiotic stress 
responsive gene 16 . In salt tolerant model plant Thellungiella, the 
expression level of cyclophilin was found to be higher 44 . Role of 
CYP20-3 (At3g62030) gene from Arabidopsis has been reported 
in light and other stress conditions 16 . Another gene from rice 
OsCyP20-2 (LOC_Os5g01270), showed tolerance against various 
abiotic stresses 16 . In the case of yeast cyclophilins, it has been also 
reported that CPR1 gene help in countering multiple abiotic stresses 
when overexpressed in E. coli and Yeast 16 . The cyclophilin protects 
the cells against oxidative induced death 42 and expression levels were 
found to be elevated after oxygen stress 43 . However, the role of cyclo- 
philin from P. indica, a salinity tolerant fungus, in stress condition 
has not been described so far. Here we have dissected the structure- 
function of PiCypA in osmotic stress tolerance in tobacco plants and 
demonstrated its novel RNA-binding property. 

The PiCypA overexpressing transgenic tobacco plants exhibited 
better growth and survival ability as compared to normal WT and 
VC plants under osmotic stress. This study suggests that PiCypA 
may play a crucial role in stress protection presumably through its 
active role in maintenance of protein structure during stress. The 
overexpression of PiCypA in tobacco plants leads to improvement in 
the antioxidant machinery and accumulation of osmoprotectant like 
proline, which cope up with the damage caused by elevated ROS 
during the stress. The harmful effects of increased ROS and 
improved antioxidant machinery in transgenic rice under stress have 
also been reported earlier 8,21,45 . PiCypA expressing T 2 transgenic 
plants also showed normal growth and increased osmotic stress tol- 
erance indicating the functional introduction of the trait, which is 
stable in T 2 generation of transgenic tobacco plants. Interestingly, we 
have found several interacting partners of PiCypA from P. indica and 
plant homologues of few of them have been reported in abiotic stress 
response previously such as PGK, HSP6 and eIF-4A 9,33 36 . The 
expression levels of these genes were significantly higher in 
PiCypA overexpressing lines as compared to normal WT plants. 

Although symbiotic relation of P. indica fungus is known to pro- 
vide salt stress to plants still molecular mechanism of such action is 
currently unknown. To understand the molecular mechanism we 
determined the atomic-resolution structure of a cyclophilin A-like 
protein from P. indica and observed novel RNA binding property. 
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Figure 4 | Expression level of tobacco homologous genes of PiCypA interacting partners (PGK; phosphoglycerate kinase, HSP6; heat shock protein 6, 
eIF-4A; DEAD box helicase). (a), (b) & (c) Real time expression analysis of PGK, HSP6 and eIF-4A, respectively (where WT-N is normal WT plant 
and WT-S is WT plant with stress treatment in NaCl). (d) Semi-quantitative PCR of these genes PGK [1], HSP6 [2], eIF-4A [3] and actin [4] as control. 
The gel is cropped from four gels and all four gels are shown in supplementary figure S5. 



Crystal structure of PiCypA revealed a conserved canonical cyclo- 
philin like fold comprising two a-helices (P33 - K44 and M139 - 
V148) and eight y3-strands (V9 - 115, Q18 - L27, H57 - 160, 
F63 - G67, G99 - M103, F115 - T118, F132 - E137 and V159 - 
T166) in (j8i-2)-<Xi-(/?3-5)-#2-/?8 architecture. In addition to usual 
cyclophilin like properties we also observed a novel RNA binding 
property in PiCypA. As CsA titration data indicates (Supplementary 
fig. S3a-c), perturbed residues upon RNA interaction are shared 
with those forming CsA interacting region (M103, A104, W124, 
K128, H129 and V130). All residues involved in RNA interaction 
are conserved in cyclophilins from different organisms throughout 
the evolution. PiCypA like AtCyp59 has specificity for GU rich 
sequences indicates a possibility of its role in stability of mRNA 
transcripts as —70% of cellular transcripts carry UG-rich region in 
their 3'-UTR as well as exon-intron junctions as discussed earlier 46 . 
Other Cylophilins like AtCyp59 from Arabidopsis thaliana and 
hCyP33 from human have been reported to bind RNA but they 
contain an additional RNA recognition motif (RRM) 46,47 . Though 
PiCypA does not have any canonical RRM domain it is still able to 
bind RNA provides an interesting insight into a new class of RNA 
binding motif. Possibly PiCypA has evolved and got rid of RRM 
domain but RNA binding property remains conserved like other 
RNA binding cyclophilins 46 . PiCypA might be involved in stabiliza- 
tion of RNAs and possibly protecting the RNA integrity by chaper- 
one-mediated process. This opens up window for further exploration 
for identification of target RNAs that are interacting with PiCypA. 

In this study we have reported the structure of a protein from P. 
indica fungus at atomic-resolution. We have also demonstrated for 
the first time a direct evidence of countering osmotic stress tolerance 
in plant by genetic modification using a P. indica gene. A novel RNA 



binding activity of PiCypA was also observed that might be part of 
machinery involved in RNA stabilization during stress conditions. 
RNA binding activity of PiCypA generates new avenues for involve- 
ment of cyclophilins in various mechanisms involving RNA stability 
and other processes involving transcripts and hnRNPs. These find- 
ings open up new areas for investigating the cellular machinery in 
plants and fungus to counteract salinity stress tolerance and have a 
great potential for crop improvement. 

Methods 

Identification and cloning of cyclophilin A-like gene of P. indica. Identification of 
candidate genes responding to osmotic stress was done as per the method previously 
described 30 . One of the clones, cyclophylin A-like protein (accession number 
GQ214003), was selected for further studies. The gene cloned into pBSK vector was 
amplified by using PCR Primers (forward: 5' CTCGAGCATATGTCCCAGCCCAA 
CGTCTACTTTG 3' and reverse: 5'-GAATTCTTAGACAGTGCCAGACGCAGT 
AATCTTG 3'). The PCR product (495 bp) was cloned into the pGEMT-easy vector 
(Promega) and sequenced using the T7 and SP6 primers respectively. Subsequent 
subcloning was done in pET28a vector (Novagen) using the Ndel and EcoRI 
restriction sites, generating the pET28a-PiCypA construct. 

Expression and purification of PiCypA. The native PiCypA was expressed in E. coli 
BL21 (DE3) codon plus cell in LB medium. For NMR samples preparation, E. coli 
BL21 (DE3) codon plus cells were grown in 15 NH 4 C1 and l3 C-Glucose labeled media. 
Cell cultures were grown at 37 C till OD 600 ~ 0.8, induced with Isopropyl |3-D- 
thiogalactoside (IPTG) overnight at 18 'C, and then harvested by centrifugation. 
Soluble fraction of lysates was transferred onto a 5 ml Ni-NTA affinity beads 
(QIAGEN) on a gravi column (GE Healthcare), pre-equilibrated with lysis and wash 
buffer. Elutes containing PiCypA were buffer exchanged and concentrated against 
50 mM NaCl, 20 mM Tris-Cl buffer pH 6.5. Purified protein was then subjected to 
Thrombin (Amersham Biosciences) digestion to remove the N-terminal His-tag. 
Traces of thrombin were removed by passing the reaction mixture through a 
Benzamidine column (GE Healthcare). Further purification of PiCypA was carried 
out by size-exclusion chromatography by HiLoad 16/60 Superdex 75 preparative 
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column (GE healthcare). Purified fraction was concentrated to ~1 mM (18 mg/ml) 
of protein concentration using 3 kDa cut-off Centricon (Millipore). 

Plant transformation. The PiCypA from pGEMT digested by Xhol and EcoRI and 
subcloned into pRTlOl and from this the complete cassette (35S CAM promoter + 
cyclophilin ORF+ NOS terminator) subcloned into pCAMBIAl301 vector. The 
pCAMBIA1301-P/CypA also carries the HYG (Hygromycin) gene as a selectable 
marker. Tobacco (Nicotiana tabacum cv. tabbacci) leaf disks were transformed 
through a procedure 9 with Agrobacterium tumefaciens (LBA4404) containing 
pCAMBIA1301-P/CypA. Putative T 0 transgenic plants were regenerated from the 
callus in the presence of Hygromycin and further screened by using PCR and |i- 
glucuronidase (GUS) assay 48 . The seeds from these plants, i.e., T 0 seeds, were 
germinated on Hygromycin containing medium to select for all of the transgenic T 1 
seedlings. GUS activity was visualized in the leaf tissue and fluorescence emission 
intensity was measured using 1 mM 4-methyl-umbelliferyl-|3-D-glucuronide as 
substrate 48 . 

Leaf disk assay. Leaf disks of 1.0-cm diameter were excised from healthy and fully 
expanded tobacco leaves of similar age from transgenic and WT plants (60 days old). 
The disks were floated in a 4 -ml solution of 200 mM NaCl and water (experimental 
control) for 72 h 49 and then used for measuring chlorophyll a and b 
spectrophotometrically after extraction in 80% acetone 50 . The experiments were 
repeated at least three times with different transgenic lines. 

Genomic DNA PCR and Southern analysis. Genomic DNA was isolated by 
grinding leaf tissue followed by extraction using the CTAB (N-acetyl-N,N,N 
trimethylammonium bromide) method 51 . PCR was performed using the gene specific 
primers PiCypA-Fl (5' CTCGAGCATATGTCCCAGCCCAACGTCTACTTTG 3 ' ) 
starting from the translation initiation site and PiCypA-Rl (5'-GAATTCTTAGAC 
AGTGCCAGACGCAGTAATCTTG 3'). For Southern analysis, 20 ug of genomic 
DNA from PCR-positive tobacco lines was digested with Hiwdlll, electrophoresed, 
and blotted on Hybond N membranes (Amersham Pharmacia). Radiolabeled 495 bp 
complete ORF of PiCypA cDNA was used as a probe. Hybridization and washing were 
carried out at 55°C 52 . 

Real time PCR, immunoblotting and in vitro RNA binding assay. Total RNA was 
extracted from 200 mg of young leaves of tobacco transgenic lines using the RNeasy 
total RNA isolation kit (Qiagen, Germany). RT-PCR was performed using total RNA 
pre-treated with RNase-free DNase L 5 ug aliquot of total RNA from each sample was 
used in a 10 ml reverse transcription reaction system with oligo-dT as primer 
(Takara, Japan). 1 ul aliquot of reverse transcription reaction mixture after the 
completion of cDNA synthesis was used as the template for real time PCR 
amplification. Total soluble proteins from the tissue were separated by 12% SDS- 
PAGE, and analysis was performed by using PiCypA specific primary antibodies 
( 1 : 5,000 dilutions) and alkaline phosphatase conjugated secondary antibodies 52 . The 
in vitro RNA binding assay was done using the method described previously 53 . For 
this 1 ug of BSA and increasing concentration of PiCypA (1, 2 and 4 ug) proteins 
were dot-blotted on precharged PVDF membrane and probed with 13 mer oc or y 32 P 
radio-labeled RNA oligo (5'-AUAGCCUCAACCG-3 r ). 

Root growth assay. Seeds of the T 0 transgenic tobacco plants were selected by 
germinating them on Murashige and Skoog 54 medium in the presence of hygromycin. 
These seedlings representing the Ti generation were tested for their sensitivity of root 
growth to osmotic stress. Individual seedlings were selected for uniformity and their 
initial root lengths were recorded. 

Biochemical and physiological parameters. Quantitation of H 2 0 2 , lipid 
peroxidation, electrolyte leakage, antioxidant enzyme assays and proline estimation, 
were performed using the method described previously 55 . Relative water content 
(RWC) was calculated following the method of Barrs and Weatherley 56 . 

Crystallization, data collection and processing. Purified PiCypA protein was used 
for crystallization. After screening and optimization diffraction quality crystals 
belonging to crystal space group C222i were obtained. The condition of 
crystallization and data collection has been reported earlier 30 . Diffraction data were 
collected with Cu K-radiation (X = 1.5418 A) using a Rigaku FR-E+ SuperBright 
microfocus rotating- anode dual -wavelength (Cu and Cr) X-ray generator equipped 
with R-AXIS IV + + detectors. The crystals diffracted to 1.97 A resolution limit. The 
data sets were indexed, integrated and scaled using MOSFLM 57 and ccp4 suite 58 . 
Phasing was obtained by molecular replacement 59 using human Cyclophilin A 
Structure (Protein Data Bank ID code 3K0N) as a search model. Model building, 
phasing and refinement were performed using AutoBuild, phenix-refme and RefMac 
software on Phenix and ccp4 suits 59 . Detailed crystallographic statistics are listed in 
Supplementary Table SI. 

NMR spectroscopy. Spectra measurements were carried out at 30 °C on Bruker 
Avance III spectrometers equipped with 5 mm cryogenic triple resonance TCI probe- 
head, operating at field strengths of 500 and 700 MHz. Standard set of triple 
resonance was used to obtain sequence-specific backbone and side-chain resonance 
assignments. Details of resonance assignments have been reported earlier 22 . 
Temperature calibration was performed using a 100% d 4 -methanol sample 60 . 
Backbone dynamics of PiCypA was studied using 15 N relaxation data (R lt R 2i and 



15 N-f H} nOe) obtained from 2D [ 15 N, 'H] NMR spectroscopy 61 " 63 . The relaxation 
data were analyzed using the MODELFREE approach 64-65 . 

All titration experiments were performed by adding NaCl, Cyclosporin A or RNA 
in equal molar steps and 2D [ 15 N, l ti] HSQC spectrum was measured at each molar 
ratio till no further chemical shift changes were observed. The cumulative backbone 
chemical shift perturbation (CSP) was calculated using the formula: 

A5( 15 N, 1 H n )=SQRt{[5A5( 1 H n )] 2 +[A8( 15 N)] 2 } 

where AS is the change in backbone amide peaks chemical- shift of PiCypA in the 2D 
[ 15 N, 'HJ-HSQC spectra upon complex formation. Topspin2.1 (Bruker AG) was used 
for processing of all spectra and data analysis was performed with CARA 67 . 

Isothermal titration calorimetry. Isothermal titration calorimetry (ITC) 
experiments were performed by injecting 1.2 mM RNA in the cell containing 
-50 uM of PiCypA on MicroCal ITC 20 o (GE Healthcare). PiCypA was in 10 mM 
NaCl, 10 mM phosphate buffer (pH 6.5) and RNA was reconstituted in nuclease free 
water. Data analysis was carried out on Origin 8 software (OriginLab). 

Yeast-2-hybrid assay. A Gal4 based two-hybrid system was used as described by the 
manufacturer (Clontech, USA; http://www.clontech.com/). The coding region of the 
PiCypA (495 bp) was amplified by PCR with primers harboring restriction sites and 
cloned in-frame into the Ndel and EcoRI sites of the binding domain vector pGBKT7. 
The resulting vector pGBKT7-PiCypA was co-transformed with P. indica AD library 
into yeast strain AH109 harboring two reporter genes (HIS3 and /?-galactosidase) by 
the lithium acetate method 68 . 
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